% Ejemplo de invocación
% rk4(0,48,10,48,"f1");
% como ultimo parametro recibe la función del modelo
% estan disponibles f1 y f2 (el que tiene u)

function energy=Puntod2()
% Parametros de Runge kutta-4
a = 0;
b = 48;
ya = 10;
M = 48;
func = "f2";

% modelo de simulación
kArr=[0.25,0.5,1,2,5,10,15,20];
global ku;
global k;
global TOUT_M;
global TD;

TD=20;
TOUT_M=10;
ku=0.7;
colors = {"1", "2", "3", "4", "5", "1*", "2*", "3*"};
labels = {'k=1/4','k=1/2','k=1','k=2','k=5','k=10','k=15','k=20'};

% Curvas de la temperatura T(t) para distintos valores de k
for i=1:length(kArr)
	k=kArr(i);
	RK = rk4gral(a,b,ya,M,func);
	energy(i,1)=k;
	energy(i,2)=uEnergy(RK);
end

end


function energy=uEnergy(RK)
	global TD;
	energy=0;

	T=RK(:,2);
	for i=1:length(T)
		if(T(i) < TD)
			energy=energy+4;
		else if(T(i) > TD+1)
			energy=energy+1;
		end
	end
end
